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We look for possible spectral features and systematic effects in Fermi-LAT publicly available high- 
energy gamma-ray data by studying photons from the Galactic centre, nearby galaxy clusters, nearby 
brightest galaxies, AGNs, unassociated sources, hydrogen clouds and Earth Limb. Apart from 
already known 130 GeV gamma-ray excesses from the first two sources, we find no new statistically 
significant signal from others. Much of our effort goes to studying Earth Limb photons. In the energy 
range 30 GeV to 200 GeV the Earth Limb gamma-ray spectrum follows power-law with spectral 
index 2.87 ± 0.04 at 95% CL, in a good agreement with the PAMELA measurement of cosmic ray 
proton spectral index between 2.82-2.85, confirming the physical origin of the Limb gamma-rays. In 
small subsets of Earth Limb data with small photon incidence angle it is possible to obtain spectral 
features at different energies, including at 130 GeV, but determination of background, thus their 
significances, has large uncertainties in those cases. We observe systematic 2a level differences in 
the Earth Limb spectra of gamma-rays with small and large incidence angles. The behaviour of 
those spectral features as well as background indicates that they are likely statistical fluctuations. 



I. INTRODUCTION 

The existence of cold dark matter (DM) of the Uni- 
verse pQ is verified beyond any reasonable doubt [2] ■ The 
leading paradigm is that the DM consist of weakly in- 
teracting massive particles (WIMPs) whose thermal relic 
abundance is predicted to be around the observed value 
if WIMP mass is 0(100) GeV. Today the nature and 
properties of the DM particles are not known. In com- 
parison with other search methods for the DM, some- 
what unexpectedly the indirect searches for DM anni- 
hilations/decays in cosmic rays have generated lots of 
activity in the field during last few years [3]. 

As a new development in DM indirect searches a 
130 GeV gamma-ray line [H [5] has been observed in 
Fermi-LAT [6] publicly available data first from the 
Galactic centre [7HTU] and then from the six nearby 
galaxy clusters [TTJ [TJ] consistently with Fermi-LAT con- 
straints [T3]. In addition, an evidence for a double peak 
structure in photons from Fermi unassociated sources has 
been claimed [T5]. However, this result has been criti- 
cized by analyzing the low energy [1IJJ [20] as well as the 
high energy [21] Fermi data. Both the Galactic centre 
and the galaxy clusters are known to be DM dominated 
objects and the most promising known places for searches 
for DM signals. The coincidence of the two excesses both 
in shape and in energy suggest that they both originate 
from DM annihilations (disfavouring DM decays due to 
large boost factor in galaxy clusters) [TT]. The expected 
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signal from nearby dwarf galaxies is too weak for de- 
tection with present Fermi-LAT statistics [T2J Q3] . The 
improved Fermi-LAT energy resolution and larger statis- 
tics would eventually allow one to tell whether the ob- 
served 130 GeV peak is a single monochromatic gamma- 
ray line [51 [5] , perhaps from internal brehmstrahlung [7] , 
two narrow lines from 77 and Zj final states [101 1151 116] 
or a narrow box- like spectrum [17] . The observed spec- 
tral feature is shown not to be consistent with broad 
photon spectra induced by conventional DM annihilation 
modes f22j to charged standard model (SM) particles or 
with broad box like spectrum predicted by DM annihila- 
tions via light intermediate states [53]. Additional con- 
straints occur from soft secondary spectrum induced by 
Z, H or other non-stable final state particles [211 ■ This 
result imposes severe constraints on building models of 
DM with those properties [25 62 . Future observations 
are able to distinguish between those scenarios defini- 
tively [63H66]. 

Although the observed 130 GeV gamma-ray excess in 
Fermi-LAT data is perfectly consistent with the DM an- 
nihilation scenario, one should worry whether this actu- 
ally is a signal of new physics or not. It is not excluded 
that this excess is just an upward fluctuations of the back- 
ground ]§T\ or a systematic detector effect. Although the 
Fermi-LAT detector is calibrated both with Monte Carlo 
and with real data 68J, and the available datasets are 
reconstructed accordingly, there is still a possibility that 
some systematic instrumental feature affects the photon 
spectrum at very high-energies. If one finds that the 
observation indeed represents a true signal, one should 
still be worried whether Fermi-LAT sees some new astro- 
physical phenomenon [69, 70] or new physics beyond the 
SM. Comparison of data obtained from different sources 
allows one to address those questions. 
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The aim of this work is twofold. First, we search for 
and study spectral features in gamma-rays coming from 
several known cosmological objects like nearby galaxies, 
dwarf galaxies, active galactic nuclei (AGN), molecular 
hydrogen clouds, compact high-velocity clouds, Earth 
albedo and Earth Limb in addition to the galactic centre 
and galaxy clusters. If the 130 GeV feature is also seen 
in data from any of those places, this would allow us to 
study and to understand the origin of the peak. Sec- 
ond, we search for possible detector dependent system- 
atic effects in the available gamma-ray data by comparing 
gamma-ray spectra from different sources. Our aim is to 
find out whether the 130 GeV excess could show some 
systematic instrumental featrures. When doing that we 
focus on the possible dependence of the gamma-ray spec- 
tra on spatial origin, on the incidence angle 9 and the 
zenith angle Z. The first choice is motivated by the fact 
that the 130 GeV excess is seen only in very particular 
places in the sky not in all data. The second choice is 
motivated by the possible presence of a 130 GeV spec- 
tral feature also in the Earth Limb data [7IJ[72]. If true, 
this would be a surprising result since this data is used 
to calibrate the Fermi LAT detector [5H] and no spec- 
tral features are seen by Fermi-LAT Collaboration in the 
early data [73] . Therefore a careful study of possible sys- 
tematic dependences of the available gamma-ray spectra 
on 9 and on Z is well motivated. 

We find no statistically significant spectral features in 
the gamma-ray spectra in any of the above mentioned 
new sources. Therefore we proceed to detailed studies 
of Earth Limb gamma-rays. No spectral features occur 
in the total Earth Limb gamma-ray data that follows a 
power-law with spectral index 2.87 ± 0.04 at 95% CL. 
This is in a perfect agreement with the PAMELA mea- 
surement of diffuse cosmic ray proton flux at Earth with 
the spectral index 2.82-2.85, depending on the considered 
energy range |74| . confirming the physical origin of the 
Limb photon flux. 

We note that the gamma-ray data from the Earth Limb 
is obtained in different Fermi observation modes. Sta- 
tistically dominant Earth Limb data is obtained during 
the normal survey mode of LAT and is featured by very 
large incidence (off-axis) angle 9 and the zenith angle 
110° < Z < 114°. In the normal survey mode this signal 
is considered as a very bright Earth atmosphere back- 
ground to cosmological observations and is cut out by 
requiring small values of the zenith angle. An order of 
magnitude smaller fraction of Earth Limb data is col- 
lected during special observation periods when the LAT 
is facing Earth, possibly observing objects close to the 
atmosphere. One can analyze those sets of data sepa- 
rately to search for possible systematic effects. We study 
an incidence angle and time dependences of the Earth 
Limb signal. We find that in some cases spectral fea- 
tures occur in the statistically limited gamma-ray spec- 
tra. However, in those cases also determination of the 
background suffers from large fluctuations. We observe 
some difference in spectral indices of large and small in- 



cidence angle photons from Earth Limb. This effect may 
be due to small statistics of the latter dataset since, if 
present, it should be a detector effect. We conclude that 
the observed trends in the statistically limited small inci- 
dence angle Earth Limb data contradict trends in much 
larger survey mode data, and are most likely just statisti- 
cal fluctuations. New dedicated Fermi-LAT observations 
of the Earth Limb could resolve this issue. 

Similarly, we do not observe any obvious systematic 
effect that could discriminate the 130 GeV peak photons 
from the background ones, in agreement with the find- 
ings of Ref. [75] . The 130 GeV gamma-ray peak occurs in 
some spatial regions and does not occur in others, show- 
ing no systematic features in any parameter. Although 
the observed excess over the power-law background can- 
not be explained with systematic instrumental effects, 
we still can speculate that its peak-like shape might be 
due to a systematic deficit at 105 GeV. In this case the 
deficit must be nontrivially correlated with an increase of 
charged cosmic ray fluxes from the directions that show 
the gamma-ray excess causing the detector effect. We 
find this possibility unlikely since the peak- like excess has 
also observed in a stacked data from many directions. 

The paper is organized as follows. In section II we 
present details of our data analyses. In section III we 
present our results. We conclude in section IV. 



II. DATA ANALYSES 

In the present analysis, we consider the public Fermi- 
LAT [6j photon event data of 210 weeks (from 4 Aug 
2008 to 5 Aug 2012) within energy region from 20 to 
300 GeV. We apply the recommended quality-filter cut 
DATA.QUAL = l/LAT_CONFIG = 1. We make use of 
the ULTRACLEAN events selection (Pass 7 Version 6), 
in order to minimise potential systematical errors. We 
also tested CLEAN and SOURCE events selections, hav- 
ing negligible effect to our results. The selection of events 
was performed using the 18 April 2012 version of Science- 
Tools v9r27pl. 

For Earth Limb analyses we use the zenith-angle re- 
gion 110° < Z < 114°: most of the Limb photons 
come from the region 112 ± 1°. Angle 9 is the inci- 
dence angle between the photon and the axis of the 
telescope also called the off-axis angle: for 9 = 0° pho- 
tons go directly to the telescope. The Earth Limb pho- 
tons are divided into two subsets. For the first subset, 
which are collected during the all-sky scanning survey, 
ABS(ROCK_ANGLE) < 52. All Earth Limb photons 
that come more-or-less directly to the telescope are from 
observations where ABS(ROCK_ANGLE) > 52°. When 
analyzing other regions than Earth Limb, we apply the 
zenith-angle cut 9 < 105° in order to avoid contamina- 
tion from the Earth photons, as recommended by the 
Fermi-LAT team. 

To avoid the effect of point sources we exclude photons 
that are within an energy-independent cut radius of each 
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TABLE I: Regions in the Galaxy featuring large number of 
20 < E 1 < 300 GeV photons without 130 GeV excess. 



Region 


I (deg 


) b (deg 


iV 7 


1 


18.0 


0.0 


946 


2 


25.0 


0.0 


943 


3 


331.0 


0.0 


1041 


4 


337.0 


0.0 


1137 


5 


346.0 


0.0 


933 



source. We used all (1873) sources from the LAT 24 
month catalog [75]. The cut radius is considered 0.2° 
1 1 3 j . In addition, we tested the radii 0.15°, 0.25° and 
0.5° resulting no significant effect on final results. To 
avoid contamination with Galactic plane, we exclude all 
photons with |6| < 5°. 

We use kernel smoothing method to fit the gamma-ray 
data. This method should be optimal for searching for 
peak-like spectral features. For technical details we refer 
the reader to Ref. [5]. However, we note that we have 
updated our analyses compared to Ref. [9] using new im- 
proved Fermi-LAT energy resolution and calibrated data 
taking into account angular dependences presented in 
Ref. [H5] . The errors to the fitted spectra as well as to the 
power-law fits are calculated with bootstrap. All errors 
in this work are at 95% CL (2a). 

III. GAMMA-RAY SPECTRA, POSSIBLE 
SYSTEM ATICS AND CORRELATIONS 

A. Search for spectral features from various 
sources and systematics 

We searched for the 130 GeV peak in gamma-rays 
from several sources: from nearby bright galaxies (ex- 
tracted from 2MASS galaxy redshift catalogue [77]), 
dwarf galaxies, bright AGNs [75], and from compact 
high- velocity clouds (79] [80] . From all these sources we 
stacked the signal as we did in the case of studying the 
signal from galaxy clusters [TJJ. Additionally, we looked 
the signal from Andromeda and Crab nebula. None of 
these sources yield statistically significant signal of the 
130 GeV gamma ray line. However, for completeness 
we note that in the case of nearby molecular clouds, re- 
cently studied by Fermi-LAT Collaboration [ST], we do 
observe a peak-like feature at 130 GeV with statistical 
significance less than 2a. Most likely this is nothing but 
a statistical fluctuation. 

So far the 130 GeV excess has been claimed from 
the Galactic centre, nearby galaxy clusters, from several 
small regions in the Galaxy listed in Ref. [5] and from 
some unassociated Fermi sources [T5J. It is tempting to 
associate the last two classes of sources with DM sub- 
haloes. However, those may also be just statistical fluc- 
tuations. At the same time, the overall integrated cosmic 
gamma-ray spectrum defined by 9 < 62°, Z < 105°, pre- 
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FIG. 1: Spectra of Fermi-LAT gamma-rays from the cosmic 
background (0 < 62°, Z < 105°) and from the Earth Limb 
(110° < Z < 114°), together with 2a error bands, as func- 
tions of photon energy. The fluxes are in arbitrary units, the 
corresponding numbers of photons are presented in the fig- 
ure. Best power-law fits to data together with their spectral 
indices are also presented. 



sented in Fig. [TJ with a blue line, does not show any excess 
at 130 GeV. It follows a power-law spectrum with spec- 
tral index 2.58 ± 0.03 and appears as a background to 
the observed 130 GeV excesses. The 2a error band is 
also presented in Fig. [TJ with a grey band. In addition, 
there are several bright region in the sky, some of them 
are listed in Table [TJ for comparison with the similar re- 
gions listed in Ref. [9 , that do not show any 130 GeV 
excess in the spectrum. Comparing results from all those 
places does not allow us to propose any logically consis- 
tent way to explain the 130 GeV peak with systematic 
detector effect. First, it is difficult to explain why the 
detector effect must give an excess. Second, we see no 
correlations of the the excess with any particular direc- 
tion, with the the incidence angle nor with the number of 
photons from the signal region. The 130 GeV line seems 
to be a real excess. 

Although explaining the 130 GeV excess over the 
power-law background due to systematic detector effects 
seems to be disfavoured, one can still assume that the 
excess is due to new source of high-energy gamma-rays 
but its peak-like shape is a detector effect. This assump- 
tion implies that the excess is broad but there is a sys- 
tematic deficit around photon energies 105 GeV that pro- 
duces peak-like spectrum. For DM scenarios this scenario 
would imply that most of DM annihilation modes that 
induce broad spectrum can explain the excess. However, 
similarly to the previous case, this scenario is also dis- 
favoured by non-observation of 105 GeV deficit in data 
from other sources. However, such an effect might be 
correlated with charged cosmic ray fluxes and requires 
dedicated studies. 
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FIG. 2: Earth Limb data collected by Fermi-LAT as a func- 
tion of time. The grey (light) bars correspond to N = 64738 
photons with large incidence angle 6 > 55° collected during 
the normal cosmic survey mode. The green (dark) bars cor- 
respond to N = 7467 small incidence angle photons collected 
when the detector faces Earth. 



B. Limb photons and possible systematics 
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FIG. 3: Spectra of the Earth Limb photons divided into 
four datasets by the incidence angle as functions of energy. 
The chosen incidence angle regions, the numbers of photons 
in those and best fit power-law spectral indices are all pre- 
sented in the figure. The shaded areas show 2a error-bands 
calculated with bootstrap technique. 



Earth atmosphere is a bright source of high-energy 
gamma-rays that is a background for cosmic observa- 
tions. Therefore in the cosmic observation mode the de- 
tector faces away from the Earth. Above E 1 > 1 GeV, 
the brightest region in the Earth atmosphere is its Limb 
in which the diffuse cosmic rays, dominated by protons, 
collide with the outer layer of the atmosphere producing 
gamma-ray flux. The Fermi-LAT sees the photons which 
are directed to the detector. Moving in the direction of 
the Earth itself the atmosphere starts to absorb all cos- 
mic rays and the Earth albedo has only a soft photon 
spectrum. Because of this production mechanism, the 
Earth Limb signal in gamma-rays is expected have ap- 
proximately the same spectral index as the diffuse proton 
flux at the Earth. The latter is measured with very high 
precision by PAMELA satellite yielding the following val- 
ues of spectral index; [74] . 

7p = 2.820 ± 0.004 for 30 GeV < E < 1.2 TeV, 
7p = 2.850 ± 0.016 for 80 GeV < E < 230 GeV,(l) 

depending on the proton energy ranges. 

In this work we consider Earth Limb to correspond 
to the zenith angles 110° < Z < 114°. In fact most of 
the Limb photons come from the region Z = 112 ± 1°. 
We plot in Fig. [j] the gamma-ray flux from Earth Limb 
as a function of the photon energy together with 2a er- 
ror band. The Earth Limb data follows a featureless 
power law with spectral index 2.87 ± 0.04. This is in a 
good agreement with Eq. ([I]). Since the two experiments, 
Fermi LAT and PAMELA, measure the same spectral 
index for the Earth Limb gamma-ray and for the pro- 
ton fluxes, respectively, we conclude that the production 



mechanism for Earth Limb gamma-rays is in good agree- 
ment with theoretical prediction. Notice also that both 
the cosmic background and Earth Limb spectra in Fig. [T] 
are featureless, both datasets contain approximately the 
same number of photons, and their spectral indices are 
clearly different. Thus the physical origin of the two spec- 
tra are different. 

During the normal Fermi-LAT cosmological survey 
mode (rock_angle < 52°) the detector is facing away from 
the Earth and Earth Limb photons arrive to the detector 
only at very large incidence angles. At the same time 
there exist Limb data collected when detector faces the 
Earth atmosphere, rock_angle > 52°. To demonstrate 
how the Limb data is collected over the time we present 
in Fig. [2] the time histogram of the collected number of 
photons divided by the incidence angle. Photons with 
9 > 55° vastly dominate the data (N = 64738) while 
the number of photons collected when LAT faces Earth, 
6 < 55°, is an order of magnitude smaller, N — 7467. 
The Earth Limb spectrum in Fig. [T] is entirely dominated 
by the large incidence angle photons. In the following we 
study whether the spectrum depends on the incidence 
angle or not. 

We divide the Earth Limb photons into four datasets 
by the incidence angle and plot the corresponding spectra 
against the photon energy in Fig. [3] Normalizations are 
arbitrary, the number of photons and their best fit power- 
law spectra are also presented. The shaded areas show 2a 
error-bands calculated with bootstrap technique as de- 
scribed in [9] . Clearly the statistically dominant datasets 
corresponding to 60° < 9 < 75° agrees well with the to- 
tal spectrum in Fig. [I] However, the spectra for small 
incidence angle show spectral features at high energies. 
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FIG. 4: Gamma-ray spectra for incidence angles 30° < 8 < 
45° for two zenit angle ranges 110° < Z < 112° and 112° < 
Z < 114° together with their best power-law fits. The latter 
are for photon energies E < 80 GeV. 



FIG. 5: Spectra of Limb photons with small incidence angles, 
8 < 55°, collected by Fermi-LAT during three equal time 
intervals. The best power-law fit to all 8 < 55° Earth Limb 
data is shown in dotted line. 



In particular, the 15° < 9 < 30° dataset shows a peak at 
150 GeV, the 45° < 9 < 60° dataset shows a deficit at 
105 GeV and the 30° < 9 < 45° dataset shows a peak at 
130 GeV. The existence of the latter peak has been used 
by several people anonymously (by our referees) and pub- 
licly (private discussions) to argue that its existence may 
indicate that the 130 GeV peak from the Galactic centre 
cannot be physical. Notice also that the spectral index of 
the best power-law fit to Earth Limb data with small 9 is 
systematically smaller than 2.87 ± 0.04 that is the spec- 
tral index for the large incidence angle Earth Limb data. 
This affects calculations of the statistical significances of 
the spectral features in Fig. [3] Clearly, the 130 GeV 
peak (red line) is statistically insignificant (~ 2.7a) when 
compared with the power-law fit to 30° < 9 < 45° data 
(spectral index 7 = 2.78+q; 1 7) but looks much more sig- 
nificant (~ 3.5ct) over the fit to 60° < 9 < 75° data with 
7 = 2.90 ±0.04. However, the errors are large because the 
datasets are small, and the fits in Fig. [3] are consistent 
with each other at 2a level. 

Because the statistics is poor for 9 < 60° and because 
there is no clear systematic feature shared by all the lines 
in Fig. |3j our first explanations to those findings is sta- 
tistical fluctuation. To study those features further we 
concentrate on the incidence angle region 30° < 9 < 45° 
and split the gamma-rays from that region into two 
datasets according to their origin from either from the 
inner or the outer half of the Earth Limb corresponding 
to 110° < Z < 112° and 112° < Z < 114°, respectively. 
The results are presented in Fig. |4j The two spectra have 
opposite behaviour at 80 GeV and different shapes of ex- 
cesses between 100-130 GeV. At larger photon energies 
the fluctuations become large due to the lack of statistics. 
Assuming that up to the photon energies E < 80 GeV the 
statistics is sufficient for determining the power-law fit to 
data, we also present those fits in the figure. The spectral 



indices of those are quite different but consistent within 
errors. Because no systematic difference in the behaviour 
of the Limb photons is expected for 110° < Z < 112° and 
for 112° < Z < 114°, those results indicate fluctuations 
due to limited statistics. 

To study the time dependence of the spectral fea- 
tures we split the Earth Limb data with small incidence 
angles, 9 < 55°, presented with green bars in Fig. [2] 
into three datasets according to the time the events are 
recorded. The time intervals are all equal but the number 
of recorder photons depends on the time interval. The 
resulting spectra are plotted in Fig. [5j In addition wc 
plot in Fig. [5] the best power-law fit to the total Earth 
Limb data with 9 < 55° . The resulting spectral index is 
7 = 2.76 ± 0.09 that should be compared with the spec- 
tral index of the fit to all Earth Limb data 7 = 2.87±0.04 
plotted in Fig. [T] The two agree within 2a errors. Such 
a systematic dependence of the power-law spectral index 
on the incidence angle is qualitatively observed by Fermi 
and we used the results presented in Ref. [SS] to calibrate 
data. However, we find that this is quantitatively just a 
small effect that cannot explain our finding. The trend 
that the small and the large incidence angle gamma-rays 
tend to follow somewhat different power-law is a possible 
systematic feature in the data. 



IV. DISCUSSION 

The 130 GeV gamma-ray excess from the Galactic cen- 
tre, that is statistically the most significant one so far, is 
visible only in photons coming to the detector with small 
incidence angles. Wc took the photons from Galaxy cen- 
tre with 9 > 55° and we do not observe any peak signal 
in this case. However, the number of photons for this 
data sample is relatively small and we cannot draw any 
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statistically meaningful conclusion. Studies of the Earth 
Limb photons in this work show the same behaviour: the 
gamma-rays coming to the detector at small off-axis an- 
gle may show possible spectral features demonstrated in 
Figs. [3][5] while the gamma-rays coming to the detector 
at large off-axis angle do not show any spectral features, 
see Figs. [T] and [3j However, the crucial difference between 
those two cases is that for the Galactic centre signal the 
< 55° data is statistically absolutely dominant while 
for the Earth Limb signal the data is statistically domi- 
nated by 9 > 55° photons. The total Earth Limb signal 
follows featureless power-law which spectral index 2.87 
that agrees very well with the PAMELA measurement of 
the diffuse proton spectral index at Earth. We believe 
that this coincidence is not an accident but a confirma- 
tion of the theoretical prediction that the Earth Limb 
photon spectrum must have the same spectral index as 
the dominant charged cosmic ray background. Therefore, 
the Earth Limb data confirms the known physics. 

Small fraction of Earth Limb data is collected when 
Fermi-LAT has observed objects close to the Earth at- 
mosphere. In this case photons arrive to the detector at 
small incidence angle. We studied this subset of data in 
order to find spectral features and systematic effects in 
the data. As seen in Figs. |3][5j indeed, there are spec- 
tral features at large energies. However, those spectral 
features are peaked at different energies and change their 
shapes when different subsamples of this data are consid- 
ered. The possible explanations to those findings include: 



1. The observed features in Figs. |3][5j are due to statis- 
tical fluctuations. The studied subsamples of data 
are statistically limited, there are just few pho- 
tons at large energies and one expects upward and 
downward fluctuations to appear. The described 
behaviour of the spectral features in Figs. |3][5] sup- 
ports this explanation. 

2. The small incidence angle data is collected when 
observing objects close to the Earth atmosphere. 
It is possible that those observations are contami- 
nated by photons coming either from those or other 
cosmological objects or by diffuse gamma-rays from 
the background. We do not observe any distinct 
bright regions in the Limb signal and disfavour this 
explanation. 

3. The spectral features could be systematic detector 



effects due to charged cosmic rays. For example, if 
the charged cosmic ray flux from some particular 
direction is higher than the usual background, the 
detector may systematically not register incoming 
photons. For example, one could argue that there 
is a deficit of gamma-rays with 105 GeV energy. 
However, this would imply that there are preferred 
directions in the diffuse proton flux at Earth that 
radiate more cosmic rays than others. We are not 
aware of those effects. 

Our favourite explanation is the first one. 

V. CONCLUSIONS 



We have searched for the 130 GeV gamma-ray peak 
in photons from several nearby cosmological objects and 
found no new candidates in addition to the known ones. 
We also searched for possible systematic effects in the 
130 GeV photons and could not identify any. The 
130 GeV gamma-ray signal seems to be a real excess over 
the power-law background. We studied gamma-rays from 
Earth Limb most extensively and found that the Earth 
Limb signal confirms the known physics - the flux from 
the Limb follows featureless power-law with the spectral 
index equal to the one measured by PAMELA for the 
diffuse proton flux. This is approximately the theoreti- 
cal prediction for the Earth Limb gamma-rays. We have 
studied in detail the possible presence of spectral fea- 
tures in subsamples of Earth Limb data featured by small 
photon incidence angles. The behaviour of those spectra 

an 
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indicates that the spectral features in our Figs, 
most likely statistical fluctuations due to small number 
of photons in those datasets. 

An independent analyses of Earth Limb data with sim- 
ilar conclusions on the 130 GeV peak as a systematic 
effect will appear simultaneously with this paper [82] . 
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